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^ Abstract 

^ We explore the self-organization dynamics of a set of entities by considering the interactions 



c/3 



X 



that affect the different subgroups conforming the whole. To this end, we employ the widespread 
example of coagulation kinetics, and characterize which interaction types lead to consensus forma- 



i-rt tion and which do not, as well as the corresponding different macroscopic patterns. The crucial 

O technical point is extending the usual one species coagulation dynamics to the two species one. 
O 

'""' This is achieved by means of introducing explicitly solvable kernels which have a clear physical 

►>. meaning. The corresponding solutions are calculated in the long time limit, in which consensus 

■^ may or may not be reached. The lack of consensus is characterized by means of scaling limits of 

o^ 

1-H the solutions. The possible applications of our results to some topics in which consensus reaching 

f— >) is fundamental, like collective animal motion and opinion spreading dynamics, are also outlined. 

o 
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I. INTRODUCTION 

The property that characterizes those systems which long range order is due to the cou- 
phng of simple interactions at their different scales is usually denoted as emergence. It has 
been studied in many different contexts, ranging from collective animal motion [1] to lin- 
guistic dynamics [2] and passing through condensed matter p] and financial systems |lj. It 
is thus a fundamental property of many important systems which encodes essential informa- 
tion about them. Emergence is commonly studied as arising directly from the interactions 
of the most fundamental units configuring the system. However, self-organization does not 
necessarily happen this way. Macroscopic order may be a consequence of interactions at 
different scales, in which the diverse subsystems play a particular contributing role. We 
will illustrate how this occurs by means of studying the dynamics of particular coagulation 
equations. While the origins and consequences of the arising of consensus by group level 
interactions go far beyond the use of coagulation equations, these constitute a reasonable 
approach to the subject, precisely due to their universality. Indeed, coagulation equations 
have proved their utility in topics as diverse as aerosols [5], polymerization [B], Ostwald 
ripening [3, [8j, galaxies and stars clustering [9], and population biology ^0\ among many 
others. They thus exemplify the broad applicability that the understanding of this sort of 
ordering processes may bring. 

A classical example of self-organizing systems is collective organism behavior, which has 
been studied using a broad range of theoretical techniques [TTHT3]. This includes the adap- 
tation of some of the classical models of statistical mechanics of spin systems [H] . Indeed, 
statistical physics models and methods have been successfully adapted to situations of a 
completely different nature such as opinion formation and spreading [151 [ISj- One of the 
classical models in this context is the voter model p!71 - [T9] . which has been studied within 
both fields of statistical mechanics and probability theory. Due to the large generality in 
which these and similar processes appear it seems appealing extending coagulation models 
to two species dynamics. On one hand, the extension is natural from a methodological 
viewpoint; on the other hand, it smoothly matches with previous approaches in the sub- 
ject of self-organization. Clustering has been previously studied in population dynamics 
models [20] including swarming systems ^21j, and coagulation equations have been used in 
both swarming [22] and opinion formation models [23] • We shall show in the following how 



two-species coagulation may constitute a sensible extension of previous approaches. 

There exist many biological examples of collective motion, including insects, birds, fishes... 
In these cases the interaction among individuals gives rise to a consensus which consists in 
the selection of the direction of motion. A specific example that has been extensively studied 
both from the experimental and theoretical point of view is the problem of locusts swarming. 
The experiment performed in [21] revealed that locusts marching on a (quasi one dimen- 
sional) ring presented a coherent collective motion for high densities; low densities were 
characterized by a random behavior of the individuals and intermediate densities showed 
coherent displacements alternating with sudden changes of direction. The models that have 
been introduced to describe this experiment assume that the organisms behave like inter- 
acting particles [211 ES] . Related interacting particle models have been used to describe the 
collective behavior of many different organisms and analyzing the mathematical properties 
of such models is currently a very active research area [HI [131 [H] • 

In this paper we will not try to describe the detailed behavior of any specific organism. 
Our goal will be to explore the mathematical consequences of models of organism clustering, 
that is an alternative way of describing the collective behavior of these systems. This ap- 
proach does not rule out the possibility of particle-like interactions among the components 
of the system, but it assumes that the tendency to follow the behavior of the neighborhoods 
is strong enough to produce collective effects at the level of the whole cluster instead of 
individual elements. As we shall see, this experiment suggests interesting two-species coagu- 
lation dynamics. We also expect that our approach may be of some relevance in the field of 
opinion formation, again as an idealized first approximation describing the dynamics of two 
opinions. We will present a simplistic model which otherwise has a number of interesting 
features. It has a number of simplifying assumptions: 1) only two directions of motion are 
considered, while in realistic cases there are infinite; 2) it is assumed that interactions take 
place at the group level and pair interactions are not taken into account; 3) the probability 
laws of interaction are very specific and simple in order to allow for mathematical tractabil- 
ity. Despite its simplicity, our model is still able to predict consensus formation in some 
cases and to give rise to counterintuitive results. Coagulation models have been also used 
to describe collective organism motion in |22j . 

We start considering a one dimensional spatial situation in which the particles move 
in clusters of £ individuals and /"*"(£, t) (/~(£, t)) represents the spatial density of clusters 



moving towards the right (left) in the time slot [t, t + dt). Clusters are modified when they 
collide with other clusters, in such a way that the probability distribution obeys the following 
equation of motion 

dtr{i,t)=J2 mk,j;mj)f^{k,t)r{j,t)-^m,i;k,j)f^im,t)r{i,t)], (1) 

m,kj 

where ^{k,j]m,i) is the collision kernel: it states the Poissonian probability rates with 
which a collision among clusters with k and j particles occurs and yields clusters with m 
and i particles. This equation has been derived under the assumption that the clusters 
interact rarely with colliding clusters (i. e., the interaction Poissonian rates are small), in 
such a way that a well stirred distribution of clusters results. Therefore, it will be assumed 
that the distribution of cluster sizes is uncorrelated. Then the Boltzmann Stosszahlansatz 
/^(£, A;,t) = /{*"(£, t)/f( A;, t) follows, where /2 and /i represent the two-clusters and one- 
cluster distribution function respectively (the subscript 1 has been omitted in Eq. ([l])). 
Indeed, note that assuming small Poissonian rates physically indicates that only in a small 
fraction of the collisions there is a successful interaction among the clusters. We additionally 
assume that the collision kernel is symmetric under refiections ^{k,j;m,i) = ^{j,k]i,m) 
as a consequence of spatial isotropy. A characteristic property of this sort of kernels is that 
all particles form a single cluster after a successful collision occurs. Note that Eq. ([I]) is 
a two-species coagulation equation. The mathematical structure of one-species coagulation 
equations with exactly solvable kernels is relatively well understood (see for instance [26]- 
[28]). and we will build our progress on the unexplored two-species coagulation based on 
these previous works. 

II. RANDOM KERNEL 

To study Eq. ([I]) we introduce the generating function 

oo 

F±(^,t) = ^/±(£,t)/, (2) 

1=0 

and we will consider a family of kernels of an apparent physical meaning. Our first choice is 

$(£, m;0,m + i) = $(£, m; m + £, 0) = 1/2, (3) 

^{k,j;i,m,) = 0, otherwise, (4) 



meaning that when a colhsion takes place the colhding clusters merge into a single one which 
chooses its direction of motion with equal probability. The generating function allows to 
transform the coagulation equations into the differential equations system 

dtF+{z,t) = ^F^{z,t)F-{z,t) - F+{z,t)F-{l,t), (5) 

dtF-{z,t) = ^F^iz,t)F-iz,t) - F-{z,t)F+il,t). (6) 

The number of clusters moving in each direction will be denoted as 

N^it) = F^{l,t), (7) 

and by employing Eq. (Is]) and Eq. ^ we find 

iV+(t)-iV-(t) = Co, (8) 

d,N' = -^N-{Co + N-), (9) 



for a constant Cq. Physically Eq. (pi) means that in each successful interaction between 
clusters a cluster is eliminated and, on average, the same number of clusters moving towards 
the right and towards the left are eliminated. Eq. ^ can be integrated to yield 

N-it) = ^^^^^^^ (10) 

y^) e^o*/2iV+(0)-iV-(0)' ^ ^ 

what reveals that N-{t) -^ and N+{t) -^ Cq if A^+(0) > iV~(0); if A^-(O) > iV+(0) then 
N^{t) — )■ and N^(t) — )■ — Cq. The number of particles moving in each direction will be 
denoted as 

M±(t) = 9,F±(l,t). (11) 

By differentiating Eqs. ([s]) and ^ we find the differential system 

dtM+ = -M-N+ - -M+N-, (12) 

dtM- = -M+N- - -M'N^, (13) 

which translates into the equations 

M+ + M- = Ci, (14) 

dtM- = In-{Ci -M-)- ]^M-N+, (15) 



where Ci is a constant. The first equation expresses the physical fact that the total number 
of particles is conserved. The second equation is a linear ordinary differential equation 
which, as we know the values of A^"*" and A^~, can be exactly integrated by means of the 
method of variation of parameters. Its solution reads 

, ._ ... _ 2CiiV-(0)^ + e^°*/^[2CgM-(0) - 2CiN-{0f + CoCiN-{0)N+{0)t] 

^^ 2[iV-(0)-e^o*/2Ar+(o)]2 ' ^^^^ 

which indicates us that M-(t) -^ and M+(t) -^ Ci if A^+(0) > A^-(O); if A^-(O) > A^+(0) 
then M^it) -^ and M~~{t) -^ Ci. It is clear that the particles interacting as prescribed 
by the equidistributing kernel (|3|-(|4]) merge into Cq separate groups travelling in the same 
direction in the long time limit. In other words, we say that consensus is reached and that 
the corresponding macroscopic pattern is an ordered one with all particles travelling in the 
same direction. Another characteristic of this type of evolution, as we have seen, is the small 
variation of the cluster sizes during the transient toward consensus. 

The degenerated case in which A^+(0) = A^^(O) behaves differently. Now the solution 
adopts the universal self-similar form 



/*(£.«) = ^exp 



^ ' m=l ^ ' 



uniformly in £, where M = (M"*" + M~)/2 is a constant of the dynamics and d = gcd{n|a^ 7^ 
0}, where F'^{z,0) = J2'^=o^n^^- This solution develops in the limit d^ d,t^ Ad/M and 
t K. L This self-similar asymptotic result can easily derived assuming F^{z,Q) = F~{z,0), 
which in turn implies F~^{z,t) = F~{z,t) for all times. In this case Eqs. rtSf-TOJ) can be 
exactly integrated, because they reduce to a single equation of the logistic type. Using this 
solution into the Cauchy formula 



rt«.^_ 1 f ^*fe«) 



where the line integration is carried out in the counterclockwise direction, one recovers the 



cluster density. Taking the long time limit on Eq. ( 18 ) one easily finds Eq. ( 17) after a formal 
calculation. A rigorous proof of this asymptotic result is possible even when we have just 
the equality A^^(O) = A^~(0), and F~^{z, 0) and F~{z, 0) differ otherwise. However the proof 
is long and technical, and so it will be reported elsewhere [29]. The physical consequence 
is on the other hand straightforward, and it reveals that consensus is not reached in this 
case. Contrarily, clusters with higher numbers of components are being generated at every 



moment in both directions, and this process continues for all times. Simultaneously, clusters 
with lower numbers of components tend to disappear as time evolves. Both characteristics 



can be read from the asymptotic self-similar form Eq. (17). With respect to the number 
of particles propagating in each direction M^(t) they tend to equilibrate each other, and 
they stay the same for all times if they are equal initially. This is a consequence of the 



dynamics given by Eqs. (12)-(13) when the number of clusters travelling in each direction is 
the same. The physical reason underlying this behavior is the cluster availability to cause 
collisions for arbitrarily long times, and this is precisely what produces an unbounded cluster 
growth. Alternatively, the dynamics of the random kernel can be understood in terms of the 
initial difference among clusters moving in the different directions Cq. As stated in Eq. rtsl) 
this difference is conserved for all times, and it becomes more visible as the annihilation of 
clusters, mediated by successful interactions, progresses. The final outcome is a set of Cq 
clusters travelling in the direction in which the majority of clusters was travelling initially. 

III. MAJORITY KERNEL 

We now consider a different kernel which is still integrable. Let us propose the following 
coagulation equations 

dj^ii,t)= J2 ^f^ik,t)rij,t)-J2f^irn,t)rii,t). (19) 

k+j=i -^ m 

In this case two colliding groups merge again into a single one, which direction of motion is 
this time biased. It is more likely that the generated group moves in the same direction as the 
most populated colliding cluster. Specifically, if a group of k particles collides with a group 
of j particles there is a probability k/{k + j) that the resulting group continues to travel 
in the same direction and a probability j /{k + j) that the direction shifts. This interaction 
rule apparently constitutes a stronger trend toward order than the equidistributing one, but 
however its consequences are rather unexpected as we will see. Using the generating function 
formalism as in the former case we find 

a,(a,F+) = F-{d-^F+) - (9,F+)Ar-, (20) 

dt{d-,F-) = F+{d,F') - {d,F')N-^, (21) 



where the number of clusters A^^(t) = F^(l,t). For the number of particles M^(t) = 
d^F^{l,t) we find 

dtM+ = dtM- = 0. (22) 

In consequence, the number of particles travelling in each direction stays constant for all 
times. In this case consensus is not reached, and the macroscopic pattern corresponds 
identically to the initial disordered state for all times. This result contrasts with the former 
one, in which case an apparently weaker interaction rule leaded to order for long times. 
Physically, this effect takes place because during most of the successful interactions between 
particles, the majority forces a change in the direction of motion of the particles moving in 
the opposite direction. However, the rare collisions in which the minority wins produces a 
change in the direction of motion of several particles. At the end, the balance between these 
two effects implies that on average the total number of particles moving in each direction 
remains constant. 

In this second case we can gain further insight into the dynamics by proposing the solution 
ansatz 

F^iz,t) = -^^^[iz-l)tl (23) 

which corresponds to the real space universal self-similar form 



Ig^ (l) E 6^,, + o (i) , (24) 

^ / m=l ^ ^ 



where Q^ is the scaling function, in clear correspondence with the asymptotic solution of 



the equidistributing dynamics with the degenerated initial condition, see Eq. (17). Instead 
of giving a rigorous proof of this fact we offer here numerical evidence of it. 

To numerically compute the self-similar solution we consider the continuum equation 

dtf^ix,t) = rif^(y,t)rix~y,t)dy-f^ix,t) r riy,t)dy. (25) 

Jo ^ Jo 

By analogy with the classical coagulation equations it is reasonable to expect that the 
behavior of the discrete equation will approach the behavior of integral continuous equation 
as i/t becomes smaller. The Laplace transform of its solution is 

/■oo 

f^iz,t)= / /±(x,t)e-^^rfx, (26) 

Jo 

which is defined for Re(z) > 0. We assume the self-similar form in Laplace space 

f^iz,t)=^-cp^izt)=^-^^iO, (27) 

8 



where ^ = zt is the self-similar variable and Re(^) > 0. This corresponds to the self-similar 
form in real space 

f^{x,t) = ^^g^['^), (28) 

where the self-similar profile can be recovered from the inverse Laplace transform 

1 ^7+ioo 



^"(C) 



2Tci 



e^^if {r])dr], 



(29) 



7— JOO 



where 7 G M is large enough (larger than the real part of all the poles of v^^(-)) and we 
have used the Bromwich integral formula pO] . We note the analogy of these formulas with 



Eqs. (17)-(18) and Eqs. (23)-(24), i. e., we are again calculating the scaling form of the 



solution in the Laplace and real space formulations. The scaling functions obey the system 






[VP+-^+(0)]9g^-, 



(30) 

(31) 



as can be seen by making the substitution (27) directly into Eq. (25). By means of the new 



substitutions W'^ = (f^ — f^{0) and W = (p — v? (0) we arrive at the system of ordinary 
differential equations: 

^H/+ = W-W^, iW^^ = W+W[. (32) 

Changing variables ^ = — e"^ and ip^ = W^ -|- 1 we arrive at the four- dimensional dynamical 
system 

(33) 
(34) 
(35) 
(36) 



H^ 


= i^~H+. 


iPt 


= H\ 


h; 


= ^^H-, 


C 


= H-, 



where Eqs. (34) and (36) are actually the definitions of H^ and H respectively. This 



system is subject to the initial conditions 



^±(_oo) = l, H^{-oo)=0, 



(37) 



and which long time behavior is 



if±(+cx)) = 0, ^^(+cx)) = l-v?^(0) 



(38) 




FIG. 1: Numerical solution Ke{H ) versus r of system Eqs. (33), (34), (35), and (36) integrated 
in the complex plane. The initial conditions are Iie{H~^) = 1 — Re('i/'^) = — lO~^cos(0i) cos(^2)) 
Re{H-) = 1 - Re(V'-) = -lO-Sin(0i) cos(e2), Iin{H+) = lm{^P+) = -lO"^ cos(0i) sin(02), 
Im(i7") = Im(V'") = -10-Sin(6'i)sin(6l2), Oi = 397r/200, 02 = vr/2. 

as can be deduced from the different changes of variables that have been performed so far. 



This dynamical system defines the self-similar solutions (27) to coagulation equation (25) 



So we devote the rest of this section to showing that these solutions actually exist. In 



the dynamical system Eqs. (33)-(34)-(35)-(36) we can identify the invariant hyperplanes 
{H^ = 0} and the invariant plane {ip~^ = ip^ , H^ = H^}. All the points belonging to 
the plane {H~^ = H" = 0} are fixed points; and all the fixed points of this system belong 
to this plane. The quantity E = ip'^ip^ — H^ — H~ is a first integral of motion, and for 
the initial conditions under consideration E = 1. This means that the trajectories under 
consideration, which depart from the plane {H^ = H~ = 0} and go back to this plane in 
the infinite time limit are distributed in the hyperbola ip^ip' = 1 precisely in this limit. 



We have numerically integrated system Eqs. (33), (34), (35), and (36); for an example 



see Fig. ^ These trajectories numerically define the self-similar profiles v'^(-)- We still have 



to check that the inverse Laplace transform formula (29) makes sense for these numerical 
solutions. It turns out that the functions (p^{-) (considered as functions of the variable r, 
which is now considered to be complex) can be analytically extended to the regions Re(r) > 
L, Re(r) < 1/L, |Im(r)| < 7r/2 + e for some real L sufficiently large and some real e > 
sufficiently small (see [29]). On the other hand, we have checked numerically that in fact 
there are no singularities for the functions <f^{-) in the whole strip |Im(r)| < 7r/2 + e, Re(r) G 
(— oo, oo). This implies that the functions v'^(-) considered as functions of ^ (considered also 



10 



as a complex variable) are analytic in the half-plane Re(^) > and therefore the function 
Q^{-) can be obtained by means of the inversion formula for the Laplace transform (29). 



This numerically proves the existence of scaling solutions to coagulation Eq. (25) obeying 



the self-similar scaling (28). 



IV. CONCLUSIONS 

To conclude, the emergent properties of physical systems are determined by the interac- 
tions among its fundamental constituents, but they could also be described and understood 
in terms of the interactions among the different subgroups conforming the whole. Consen- 
sus is the microscopic characteristic underlying macroscopic order in these systems, and it 
is reached or not depending on the interaction nature. In the specific case of the major- 



ity kernel considered in this paper, Eq. (19), the absence of consensus is due to the finely 



tuned balance between the transition probabilities of the kernel and the number of elements 



changing their state in each interaction. Kernel (19) can be easily generalized to models 
with the form: 

dj^{l,t)= J2 (^Y f^ik,t)r{j,t)-J2f^ira,t)r{l,t), a>0. (39) 

k+j=l V -^Z m 

For this class of equations the finely tuned balance that takes place for a = 1 does not take 
place in general. An interesting mathematical question would be determining the range of 
values of a yielding consensus formation for these models. Notice that Eqs. (|3| and (111) 
constitute an example of interaction kernel which does not preserve the balances among 
collisions, and thus leads to ordered long-range states, despite the seemingly weaker trend 



toward order implied by it when compared to the one in Eq. (39). 

In the simple cases we have studied, consensus takes the form of one or several clusters 
travelling in the same direction. As two directions are possible, and the system is symmetric 
with respect to them, we note the similarity of this process with those systems with two 
symmetric absorbing states [221 E3] • It is clear that consensus is an absorbing state in our 
model, as collisions stop when consensus is reached. Complementarily, the active state is 
the one in which collisions still take place, which necessarily implies the existence of clusters 
travelling in both directions and the concomitant absence of consensus. We additionally note 
that for the dynamics given by the majority kernel described in the last section, consensus 

11 



is not reached in the coagulation equation approximation. So this disordered state could 
be thought of as a long lived, deterministically stable, metastable state which will finally 
decay to one of the absorbing states due to the effects of fluctuations. We have of course not 
studied such effects in the present paper, so we simply mention the possible similarity of this 
process with the dynamics of certain systems characterized by the proximity to absorbing 
states. Establishing a rigorous connection is far beyond the scope of this work. 

In the light of our results it seems possible that some self-organizing systems interact 
according to some rule which does not preserve the balance among the collisions or interac- 
tions at their different scales. Otherwise, as we have seen, if some moment is preserved then 
consensus is consequently not reached. As we have mentioned, these techniques may be ap- 
plied or adapted to swarming and opinion formation systems. Our approach constitutes an 
idealized flrst approximation to these sorts of problems. More realistic approaches should 
deal with the effect of fluctuations and flnite sizes and this way estimate how stochastic 
forces may influence ordering/disordering processes. Also, we have not taken into account 
the effect of saturation and fragmentation of clusters, which would keep the cluster masses 
finite. We have focused for simplicity on the short time behavior during which the initial 
growth of cluster masses could be described with the probabilistic rules we have considered 
herein. We leave for future work the implementation of these improvements. 
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